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ABSTRACT 

We propose an anisotropic generalisation of the line correlation function (ALCF) to 
separate and quantify phase information in the large-scale structure of galaxies. The 
line correlation function probes the strictly non-linear regime of structure formation 
and since phase information drops out of the power spectrum, the line correlation func¬ 
tion provides a complementary tool to commonly used techniques based on two-point 
statistics. Furthermore, it is independent of linear bias as well as the Gaussian vari¬ 
ance on the modulus of the density field and thus may also prove to be advantageous 
compared to the bispectrum or similar higher-order statistics for certain cases. For 
future applications it is vital, though, to be able to account for observational effects 
that cause anisotropies in the distribution of galaxies. Based on a number of numerical 
studies, we find that our ALCF is well suited to accomplish this task and we demon¬ 
strate how the Alcock-Paczyiiski effect and kinematical redshift-space distortions can 
in principle be measured via the ALCF. 

Key words: Line Correlation Function, Phase Statistics, Bispectrum, Redshift-Space 
Distortions, Large-Scale Structure 


1 INTRODUCTION 

Galaxy surveys map a growing fraction of our Universe at 
ever increasing accuracy and are routinely used to test cos¬ 
mological models. So far, agreement with the ACDM model, 
a spatially flat universe primarily filled with dark matter and 
dark energy in the form of a cosmological constant, is ex¬ 
cellent. Among their key science questions are the nature 
of dark matter, dark energy, and the concrete model of the 
early Universe during which primordial density fluctuations 
were generated. A range of next generation experiments, 
such as DESI (Flaugher & Bebek 2014), Euclid (Laureijs 
et al. 2011) and SKA (Dewdney et al. 2009), etc. are de¬ 
signed to shed light onto these questions. 

To make progress, it is crucial to optimize the amount 
of information extracted from the raw data that these exper¬ 
iments produce. Up until recently, the majority of the litera¬ 
ture focused on two-point statistics, i.e. the power spectrum 
P or the two-point function C to analyse the observed dis¬ 
tribution of galaxies. If the density field were Gaussian, two- 
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point statistics would be sufficient to specify the statistical 
properties of matter fluctuations. However, during the non¬ 
linear stage of gravitational interaction, where structures 
like filaments, clusters and eventually galaxies form, any ini¬ 
tial density field becomes highly non-Gaussian. Hence, by 
focussing on two-point statistics, a considerable part of the 
available information is ignored. 

Consequently, new statistical measures have received 
significant attention. Besides computing higher-order cor¬ 
relations (e.g. Peebles & Groth 1975; Peebles 1980; Frieman 
& Gaztanaga 1994; Scoccimarro et al. 1998), it was sug¬ 
gested to employ genus statistics (Gott, Dickinson & Melott 
1986; Hoyle et al. 2002; Hikage et al. 2002), Minkowski 
functionals (Mecke, Buchert & Wagner 1994; Hikage et al. 
2003) and others. However, all of these measures suffer from 
the conceptual limitation that they are correlated with the 
two-point function. Thus, an intriguing idea is to exclu¬ 
sively analyse information not already contained in two- 
point statistics, i.e. phase information. 

Correlations in the phases of Fourier coefficients emerge 
as a consequence of non-Gaussianity; since two-point statis¬ 
tics depend on the amplitudes of Fourier coefficients only, 
they are blind to phase information. The phase factors of 
the matter density field have been studied many years ago, 
but much of the original work in Ryden & Gramann (1991); 
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Soda & Suto (1992); Jain & Bertschinger (1998) focused on 
the evolution of single phases away from their initial val¬ 
ues. A first statistic that described the phase difference be¬ 
tween neighbouring Fourier modes was presented in Coles 
& Chiang (2000); Chiang (2001); Watts, Coles & Melott 
(2003). Shortly thereafter, their considerations were gen¬ 
eralised by Matsubara (2003), who quantified phase infor¬ 
mation by computing the joint probability density of phase 
factors, which illuminated its relation to higher-order spec¬ 
tra (see also Hikage, Matsubara & Suto 2004; Matsubara 
2006). In Hikage et al. (2005), the probability densities were 
applied to a galaxy catalogue from the SDSS DR-2, giv¬ 
ing constraints on bias models. Recently, a new measure for 
phase information, the line correlation function, was intro¬ 
duced by Obreschkow et al. (2013) as a spherically-averaged 
three-point correlation of phase factors. In the follow-up 
paper (Wolstenhulme, Bonvin & Obreschkow 2014) it was 
computed perturbatively and related to the bispectrum and 
power spectrum at leading order. 

The aim of this paper is to extend the line correlation 
function to account for line-of-sight and transverse distance 
scales separately. The introduction of such an anisotropic 
line correlation function opens up the possibility to detect 
anisotropies and distortions along the line-of-sight, which 
are inherently present in any galaxy survey due to redshift- 
space distortions. 

The paper is organised as follows: after a brief review of 
the definition and properties of the line correlation function 
in Sec. 2, we present and test a modihcation in Sec. 3, which 
is capable of quantifying anisotropies. In Sec. 4 this new for¬ 
malism is applied to simple mock fields and we investigate 
the sensitivity to the Alcock-Paczyhski effect and kinemati- 
cal redshift-space distortions in comparison to results based 
on the two-point function. We conclude in Sec. 5. 


2 THE ISOTROPIC LINE CORRELATION 
FUNCTION 

In this section, we review the isotropic line correlation func¬ 
tion (ILCF), which serves as the starting point for our gener¬ 
alisation to the anisotropic line correlation function (ALCF) 
in Sec. 3. 


2.1 Definition 

Given a density held J(a:), the ILCF measures correlations 
in its Fourier phase factors which are obtained by dividing 
amplitudes in Fourier space^, 

[■/ \ FT „ whitening 5 k IFT . , 

d(a:) —>■ Sk —Cfe = y— —>■ e{x) . (1) 

|Ofe| 

The whitened density held e(x) is thus devoid of any in¬ 
formation probed by two-point statistics. Consequently, the 
simplest measure of phase information must be based on 
the three-point correlator and since the whitening transfor¬ 
mation tends to collapse elongated objects to thin line seg¬ 
ments, as was observed in Obreschkow et al. (2013), these 
three points are chosen to be distributed equidistantly on a 


^ The zero mode is conventionally set to zero, i.e. ek=o = 0- 


straight line. More precisely, and expressed in Fourier space, 
the ILCF is dehned by 


l(r) = 




D\ 3/2 


V 


HI 


d^fci d^k2 d^fca 


(27r)3 

|felMfc2l.|fe3K2'ir/r 


where D denotes the dimension of space, V the volume of the 
survey and {) an ensemble average over all realisations of the 
density held. In writing Eq. (2) we adopted the same conven¬ 
tions as employed in Wolstenhulme, Bonvin & Obreschkow 
(2014). 

The ILCF differs from the conventional three-point cor¬ 
relator of whitened density helds in two respects. Firstly, it 
is multiplied by a prefactor (r^ and secondly, each 

held is convolved with a top-hat window function that cuts 
off high-frequency modes |fc| > 2'k jr. These corrections are 
introduced to regularize the integral in the regimes of large 
and small scales, respectively. As discussed in Obreschkow 
et al. (2013), Fourier modes which are associated with scales 
larger than the longest physical correlation length or scales 
smaller than the lowest characteristic length obtain random 
phase factors that would cause the ILCF to be dominated 
by random noise. 

It is important to reiterate that the ILCF is indepen¬ 
dent of amplitude information and for that reason not just 
a regularized integral of the bispectrum. It depends on the 
quantity {ek^ek^ekA = {exp[i{ekj^ + Ok^ + OkA])^ which 
can be determined by using the probability density function 
(PDF) of Fourier phases, P{0}. This PDF was computed 
perturbatively for mildly non-Gaussian fields by Matsubara 
(2003), showing that it is not solely related to the bispec¬ 
trum but a progression of all higher order spectra such as the 
trispectrum etc. Taking this result, Wolstenhulme, Bonvin & 
Obreschkow (2014) showed that the ILGF can be expressed 
to lowest order as follows: 


l{r) 



II 


3D/2 

1^1 Mfc2Mfcl+^2 l^27r/r 


d^fci d®fe2 


X 


B{ki,k2) 

VD(|fcT|)>(|fe2|)p(Tfei + fezi) 


cos [(fcl 


fcz) • r ], (3) 


where P(|fc|) and B{ki, k 2 ) are the ordinary power spectrum 
and bispectrum. Accordingly, we can think of the ILCF as 
compressing information from all higher order spectra into 
a single function. Only in the case of weakly non-Gaussian 
fields does the dominant contribution stem from the bispec¬ 
trum. 

For evaluating the ILCF on a given density held, Eq. (3) 
is not particularly useful. However, we arrive at a simple 
prescription to extract the line correlation by assuming er- 
godicity and replacing the ensemble average in Eq. (2) by an 
average over all translations and rotations. Integrating out 
the translation vector x as well as the orientation of r thus 
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gives 


l{r) 





\ki\,\k2\,\ki + k 2 \^ 27 T/r 


d°k2 


X Wodfei — k2\r) 


^fc2 ^—ki—h2 


( 4 ) 


l^fcl ^fc 2 fcl—#22! 
where we have defined the rotational average of the expo¬ 
nential Fourier factor as 

( cos (fcr), if D = 1 , 

= 1 Jo(fcr), if 75 = 2, (5) 

1 sin (fc r)/(fc r), if D = 3 . 


WD(kr) = (e' 


A discretized version of Eq. (4) will be used for all numerical 
computations throughout this paper. Details on the imple¬ 
mentation are given in App. A. 


2.2 Properties and Applications 

The morphology of the cosmic web, composed of structures 
such as clusters, filaments and voids, is largely encoded in 
the phases of Fourier coefficients (for a demonstration see 
e.g. Chiang & Coles 2000). Against this background, it 
is natural to assume that the ILCF is sensitive to these 
kinds of objects and, owing to its linear configuration, espe¬ 
cially to filaments. In fact, after a more systematic analysis, 
Obreschkow et al. (2013) conclude that the ILCF is gener¬ 
ally a measure for aspherical structures on scales ~ 2r. This 
holds true for oblate and prolate objects alike, so that l{r) 
is able to probe both cosmic sheets and filaments. 

Furthermore, they find two scaling relations depending 
on the size and density of the substructure. Letting ro de¬ 
note its characteristic scale and N the number of objects^, 
their filling factor is given by / = NVo/V where Vo is the 
volume of a D-dimensional sphere of radius ro- A change 
of the characteristic scale or the number of objects has the 
following effect on the ILCF: 

а) \ ^ i(r)^l(r/a), (6) 

/ = const. J 

б) N^aN ^ (7) 

Va 

The first relation specifically targets the behaviour of l{r) 
under a rescaling of the substructure and therefore should 
not be confused with a scaling law resulting from power law 
clustering. This can be derived from Eq. 3 and assuming the 
hierarchical model (e.g. Bernardeau et al. 2002), so that we 
have D ~ P^. For a power law spectrum P{k) ~ fc" we then 
get 

D+n 

l{ar) = a 2 ^ (g) 

implying that the slope of l{r) is less steep than that of the 
two-point function, which would scale as for the 

same power spectrum. This conclusion was already found 
from numerical observations in Obreschkow et al. (2013). 
Furthermore, compared to conventional correlation func¬ 
tions, the second relation is also somewhat unintuitive: in¬ 
creasing the number of objects (and hence the clustering) 

^ All Objects are assumed to have random positions and orien¬ 
tations. 


causes a decrease in the amplitude of l(r) instead of an in¬ 
crease, as one would expect. The cause of this feature is 
that the increased number of different translations and rota¬ 
tions of objects randomize the phase factors, hence decreas¬ 
ing their correlation. This effect can also be understood by 
employing the halo model for a particular class of objects 
with orientation </>, mass m and density profile ShC \ 
with Fourier transform For random positions of 

the halos, one only needs to consider the one-halo term, so 
that, according to Smith & Watts (2005), the power- and 
bispectrum become 


P{\k\) 

{2tv)o 

B(ki,k2, fes) 

{2tv)o 


J J Ql 


\Uk{4>,rn)f, (9) 




i=l 


( 10 ) 


where p denotes the total mass density of all halos and Qd 
the surface area of the D-dimensional unit sphere. If all halos 
have identical mass M, the mass function is given by n(m) = 
N/V S^{m — M) and we see that P oc 1/N, i.e. B oc l/A^. 
Hence, the combination of power spectra and the bispectrum 
in Eq. (3) gives rise to the scaling ~ 1/y/N. 

Being a relatively new tool in large-scale structure 
statistics, the ILCF has been applied only a few times so 
far. For instance, in Obreschkow et al. (2013) it was shown 
to be a significantly more sensitive tool to analyse proper¬ 
ties of dark matter than the two-point function. The au¬ 
thors examined a series of cosmological N-body simulations 
with different warm dark matter masses, finding that l{r) 
distinguishes those masses out to scales approximately five 
times larger than the two-point function. This feature is a 
direct consequence of the ILCF’s sensitivity to elongated 
structures and its ability to exclusively probe the non-linear 
stages of gravitational evolution. In another recent applica¬ 
tion (Alpaslan et al. 2014), the ILCF was used to identify a 
possibly new class of structure in the matter distribution, de¬ 
noted as tendril galaxies. These galaxies were spotted in the 
Galaxy and Mass Assembly (GAMA, Driver et al. (2009)) 
catalogue as remnants after all filament and isolated galax¬ 
ies had been removed by appropriate detection algorithms. 
They appear in thin chains that extend into voids and dis¬ 
play a different line correlation than those galaxies residing 
in the larger filaments. Finally, in Wolstenhulme, Bonvin & 
Obreschkow (2014), the ILCF was computed using tree-level 
perturbation theory giving good agreement with direct es¬ 
timates from N-body simulations. An analytic formula for 
l{r) allows it to be modelled properly and opens the door 
for measuring various cosmological parameters. 


3 THE ANISOTROPIC LINE CORRELATION 
FUNCTION 

Our ability to measure cosmological parameters with the line 
correlation function is intimately linked to our understand¬ 
ing of how it is influenced by effects that cause anisotropies 
in the distribution of galaxies (see Sec. 4). Not knowing the 
systematics of these effects increases uncertainties in the re¬ 
sult. However, an anisotropic signal does not only add com¬ 
plications, its detection would also open up new avenues to 
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Figure 1. Examples of two-dimensional density fields used 
throughout the paper. Upper panel: 50 aspherical Gaussian halos 
with (7x = 0.006 L, ay = 0.024 L and orientations drawn from a 
flat distribution. Lower panel: Zel’dovich approximation, cosmo¬ 
logical parameters are to be found in the text. 



probe and constrain the cosmological evolution in our uni¬ 
verse. For this reason, it is a vital step to extend the ILCF to 
account for radial and transverse distance scales separately, 
allowing us to detect distortions along the line-of-sight. 

Due to the corrections regularising the ILCF, there are 
various possibilities to break rotational symmetry. In the 
following we will therefore explore two different mode cut¬ 
offs, one being spherical as in the definition in Eq. (2), the 
other one aspherical. We show that the latter is better suited 
to pick up anisotropies and test some of its properties. 


3.1 A Spherical Cut-Off 

Naively, it is straightforward to write down a two- 
dimensional analogue of the ILCF, especially when adopting 
the distant-observer approximation, so that we can identify 
the line-of-sight direction with the 2 :-direction - we simply 
restrict the angular average in our prescription in Eq. (4) to 
orientations of the transverse part of the vector r, giving 

(V) // 

X cos [(fc||,i - fc||, 2 ) rii] - fex.2| rx) 

^ fill 

|5fel 5k2 S-ki-k2 I 



Figure 2. Anisotropic line correlation function with spherical 
cut-off for statistically isotropic (top) and anisotropic fields with 
von Mises parameters k = 3, fi = 0 (bottom). Each field contains 
50 halos and the width of their Gaussian profiles are chosen as 
cTx = 0.006 L and ay = 0.024 L. Results are averaged over 50 
different realisations. Colour scales are adjusted. 


Here, r± = \Jx^ and ry = z (D = denote transverse 
and radial separations, respectively. 

To see how well the function above distinguishes be¬ 
tween statistically isotropic and anisotropic density fields, 
we consider a simple test case. We construct fields consist¬ 
ing of superpositions of aspherical halos with a Gaussian 
density profile (see upper panel in Fig. 1). While their cen¬ 
tre positions are randomly distributed, their orientations are 
either drawn from a flat or a von Mises probability distribu¬ 
tion. The von Mises distribution (Mardia & Jupp 1999) is a 
close approximation of the wrapped normal distribution - the 
analogue of the usual normal distribution on a circle - and 
depends on two parameters: p, the centre of the distribution; 
and K, which is a measure for its ‘peakedness’, i.e. increasing 
K causes a higher concentration about the centre. It reads 

K cos 

p{(f)\p,K)= , (12) 

27r jq 

where Io{k) denotes the modified Bessel function of zeroth 
order. Furthermore, all fields here and in the remainder of 
this paper, if not stated otherwise, are set up in a two- 
dimensional box of sidelength L with Af = 512 grid cells 
(details are given in the captions of the corresponding fig¬ 
ures). Extracting the line correlations of these sample fields 
according to Eq. (11) leads to the plot in Fig. 2. While 
f(rx,r||) only depends on r = -|- for the statistically 

isotropic fields, as expected, we do not observe a clear de¬ 
viation from rotational symmetry in case of the anisotropic 
ones. This implies that the anisotropic signal is smeared over 
the whole (rx,r||)-plane, which renders the spherical cut-off 
either tricky or ill-suited for the purposes given in the in¬ 
troduction of this section. For that reason we consider an 
alternative method that incorporates an aspherical cut-off. 
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X 
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Figure 3. Cut-off ellipses for different pairs (rx,r||) and the re¬ 
sulting whitened and filtered density fields. Note that the blurring 
increases from top right to bottom left because of a decreasing 
number of modes. 


3.2 An Aspherical Cut-OfF 

When exploiting our freedom to alter the mode cut-off, we 
must still ensure that it regularises the integral for large 
scales and at the same time does not introduce any artificial 
anisotropies. Therefore, we impose that it should satisfy the 
following two conditions: 

(1) The number of enclosed modes scales with 

(2) The number of enclosed modes for each pair of values 
(r_L,r||) with the same r is constant. 

There are still many possible cut-offs that satisfy these con¬ 
ditions; the one that we propose here is based on oblate 
spheroids of constant eccentricity (i.e. ellipses in the two- 
dimensional (fcy, fe||)-plane). The semi-minor axes is sup¬ 
posed to scale as 27r/r and more importantly, its orientation 
is chosen to be aligned with the vector r = (ry,r||), such 
that the spheroid rotates for varying scales. This choice is 
motivated by the fact that by filtering out all modes ex¬ 
cept those that predominantly belong to a given direction 
in fc-space, objects are singled out that are aligned with the 
transverse direction in real space. This is demonstrated for 
a two-dimensional example in Fig. 3 which shows three dif¬ 
ferent cut-offs and the corresponding whitened and filtered 
density fields. The green ellipse, for instance, belongs to a 
configuration where r± = 0. Mathematically, this rotating 
cut-off may be expressed as 

0Tj{k, r) = C + (jf — l) (fe ■ r)^ , (13) 

where rj = const. ^ 1. In two dimensions, we can convert Eq. 
(13) into polar form with ifi the polar angle in the (fcy, fc||)- 



Figure 4. Ratio of variance to squared average as a function 
of radial scale r based on 20 realisations. Upper panel: density 
fields consisting of Gaussian halos with random positions and 
orientations. Lower panel: fields set up according to the Zel’dovich 
approximation, also based on 20 realisations. 


plane and we get 


k < 


2 7 r/(r 7 r) 




(14) 


-f cos ip -^ sin if 


Thus, both semi-axes scale with r~^ while the eccentricity 
of the ellipses is constant and given by e = -^/l — as 
desired. For 77 = 1 we recover the original spherical trunca¬ 
tion, while different values of rj influence the signal-to-noise 
ratio and the sensitivity to anisotropies, as we will see in the 
following. 

Since the eccentricity is left as a free parameter, we 
may ask whether there is a particular choice of rj that is 
optimal in the sense that it gives optimal signal-to-noise ra¬ 
tios for a range of scales. To determine if there is such a 
choice we employ the same sample of statistically isotropic 
fields as above, consisting of Gaussian halos with random 
positions and orientations. After computing the average line 
correlation and its variance from a number of realisations, 
we integrate out the angular dependency by averaging over 
all data points belonging to one of 25 radial bins of width 
Ar ~ 0.0016 L that extend outwards from the origin r = 0. 
Normalizing the variance by the square of the average, we 
obtain an r-dependent noise-to-signal measure that we can 
compare for different values of rj. The results of this analysis 
are shown in the upper panel of Fig. 4. All curves display a 
nearly continuous increase in the noise level, which is a con- 
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Figure 5. Same as Fig. 2 but using the anisotropic line correla¬ 
tion function with aspherical, rotating cut-off. Colour scales are 
adjusted. 
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sequence of a decreasing number of modes. However, the line 
correlation with spherical cut-off exhibits the largest noise- 
to-signal for most scales, indicating that a different choice 
for rj is preferable. Indeed, we observe that erf /Z drops when 
the cut-off becomes more elliptical and eventually reaches a 
minimum that roughly lies in the interval rj^ ~ [16,32]. For 
larger values, the noise level begins to rise again, such that 
the optimal rj must reside in the given interval. However, 
since the geometry of the cut-off matches that of the Gaus¬ 
sian halos, one might suspect that this choice of parameter 
is optimal for this particular class of density fields only. To 
check this suspicion, we employ a second, more realistic set 
of density fields which are set up according to the Zel’dovich 
approximation (Zel’dovich 1970) for a ACDM universe with 
parameters llm = 0.314, h = 0.674, erg = 0.9 and baryon 
fraction fb — 0.038 in a box of sidelength L = 1 h~^ Gpe (for 
an example see lower panel of Fig. 1). Applying the same 
analysis to the Gaussian halo fields leads to the results in 
the lower panel of Fig. 4. As before, we note that a mini¬ 
mal noise-to-signal ratio is obtained in the parameter range 
~ [16,32], strengthening the previous result^. Therefore, we 
define the anisotropic line correlation function (ALCF) for 
the remainder of this paper by adopting the rotating cut-off 
with rf — 32. After fixing rj, we still have the freedom to 
adjust the profile of the cut-off, for instance to a Gaussian 
instead of a top-hat profile, and thus minimize the noise-to- 
signal ratio even further. However, for simplicity, we work 
with the top-hat filter in this article. 

Having settled on a definition of the ALGF, let us see 
how this choice influences the sensitivity to anisotropies. To 
that end, we repeat the test case from above where the fields 


® This result was obtained in a 2D analysis; we caution the reader 
that the optimal value for r) in 3Z) might be different. 



<z> 


Figure 6. Angular dependence of the ALCF for von Mises 
dsitrlbuted density fields. Data points stem from the extraction 
procedure described in the text with error bars representing the 
l-(7 uncertainty from averaging over 100 realisations. Solid lines 
are the corresponding fitting results. 


Table 1. Input and fit parameters for von Mises distributed den¬ 
sity fields. Parameters a and /3 are fitted simultaneously to all 


data sets. 

, giving a. 

= 0.58 ± 0.02 and /3 = 

-0.23 ±0.02. 


Ki 

f^2 

K3 

A£4 

K5 

Input 

Fit 

0.5 

0.5 ± 0.1 

1.0 

1.1 ± 0.1 

1.5 

1.7 ±0.1 

2.0 

2.3 ±0.2 

3.0 

3.4 ±0.2 


are either statistically isotropic or anisotropic by giving the 
orientation of the Gaussian halos a strong preference in the 
line-of-sight direction. The resulting ALCFs are displayed in 
Fig. 5. While the isotropic fields still give rise to rotationally 
invariant line correlations, we now see a clear anisotropic 
signal in the lower panel. Since the isocontour lines appear 
squashed and the signal peaks close to ry = 0, the whole 
transverse direction is enhanced. Given that most of the 
halos point in the line-of-sight direction, this feature may 
appear counter-intuitive, but it is simply a continuation of 
the properties of the ILCF mentioned earlier. In Sec. 2.2 it 
was found that the amplitude of the ILCF decreases when 
the number of objects increases, which is reproduced here 
as well. In the following section we study the properties of 
the ALCF more closely. 


3.3 Properties and Examples 

To characterise the properties of the ALCF and verify that 
it extends those of the ILCF in a consistent manner, we 
work again with the simple Gaussian halo fields and con¬ 
sider two anisotropic setups. For the first, we introduce a 
preferred direction, as already seen in the examples of the 
last two sections; for the second, we vary the scale of the 
halos depending on their orientation. 
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(a) N« = 60 , 7Vj_ = 0 


(b) Af|| = 59, = 1 


(c) N« = 30 , Wx = 30 
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Figure 7. Line correlations for fields with different numbers of elongated Gaussians pointing either along the line-of-sight (Afy) or 
transverse direction (A^x), based on 20 realisations each. Note the change of colour scale in the middle panel. 


3.3.1 Varying Number Densities 

As before, the preferred direction is established by drawing 
the orientation angles of halos from a von Mises distribu¬ 
tion. We want to analyse how the strength of anisotropy, 
parametrised by the von Mises parameter k, affects l(r^_,r \\) 
and in particular, whether there is a scaling relation similar 
to the one valid for the ILCF. Therefore, we firstly introduce 
a probability density of objects, pxr((^), which is defined on 
the interval [—7r/2,7r/2], such that 

N X Ph{4‘) d<?!> = number of objects in the bin d<j) , (15) 

r-w/2 

^ / pjf(0)d0=l (16) 

i-7r/2 

with N, the total number of objects. Due to the symmet¬ 
ric geometry of the halos, pn is related to the probability 
density of orientation angles by 

, ,, , ,, ,, , cosh [k cos (fli — u)l 

PHW=pi^)+p{ct + 7T)^ - ^ ^ (17) 

TT lo(ei) 

where we have substituted in the von Mises distribution from 
Eq. (12). 

As already suggested by Fig. 5 and in agreement with 
the interpretation of line correlations, we expect that the 
ALCF scales inversely with the probability density of ob¬ 
jects. To verify this assertion, we compute the ALCF for 
five ensembles, each with p = 0 and a different value for k. 
Afterwards, we extract the dependence on the angle (p from 
the two-dimensional function Z(rx,r||) by averaging over a 
fixed number of points with constant (p, but constrained to 
the radial interval 0 ^ r R — 0.02 L to exclude areas of 
large sample variance, i.e. 

1 

l{(p) = — / dr Z(rsin((), rcos^). (18) 

R Jo 

The results are displayed in Fig. 6; error bars originate from 
averaging over 100 realisations. We then fit this data to de¬ 
termine the scaling behaviour and to see whether we can 
recover the respective values of k that were used as input 
parameters for the construction of the fields with the model 

l{(p\K,a,l3) = a[pH{<p\K)]^ , (19) 


where we allow for a varying amplitude and exponent. Since 
a and /3 should be universal (independent of «), we fit them 
simultaneously to all five data sets while we leave k free to 
vary amongst them. The outcome of this fit along with the 
input parameters is summarised in Tab. 1, and shown for 
four example data sets as the thick solid lines in Fig. 6. We 
note that the fitted curves reproduce well the behaviour of 
the extracted data and we recover all peakedness parameters 
nearly within the error bounds. Moreover, taking the fitted 
values of a and pi to make a prediction for fields with a 
flat probability density (i.e. pn = 1/rr), we compute I = 
0.76T0.03 while getting I = 0.753±0.001 numerically. Hence, 
we can confirm that the ALCF scales with the probability 
distribution of objects where it is likely that the parameter 
pi depends on the eccentricity of the cut-off ellipses. 

Before moving on to the second test case we briefly de¬ 
scribe another interesting feature related to the probabil¬ 
ity density of structures. To that end, we generate fields 
that can be considered as an extreme case of the von 
Mises fields above: varying numbers of elongated Gaussians 
pointing either in radial or transverse direction with fixed 
N = N±_ + A’li. From Fig. 7 we observe that a strong signal 
in the line-of-sight is received when all objects are point¬ 
ing in that direction, see panel (a). Even so, this signal can 
be altered significantly if a single object is aligned with the 
transverse direction, see panel (b). Adjusting N±_ and A^n 
until equality leads to a symmetric signal, see panel (c). 
This feature of the ALCF can be explained along the same 
lines as before: in (a), no information on structures in the 
transverse direction is present, while in (b) one object is 
present, which henceforth dominates due to minimal phase 
noise. Thus, care must be taken when using the plots of the 
ALCF to eyeball the strength of anisotropy. 


3.3.2 Angularly Varying Scales 

After having demonstrated the dependence on angular num¬ 
ber densities, we turn our attention to angularly varying 
scales. Our sample halo fields are created as follows: for each 
orientation angle (p £ [—7r/2,7r/2], we determine the scale Oy 
of the Gaussian profile from the function 

5((() I (7, Act, n) = (T + Act cos (w (()) (20) 
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Figure 8. Isocontours with l-cr error bars for two different values 
of I extracted from the ALCF of fields with angularly varying 
scales, based on 100 realisations each. Solid lines stem from the 
results of the fitting procedure. 


Table 2. Expected and fitted parameters for density fields with 
scales of explicit angular dependence. 



Rh [10-2 L] 

Rh [lo-^^L] 

A 

7 

Input 

1.14 ±0.02 

1.03 ±0.02 

0.5 

_ 

Fit 

1.16 ±0.04 

1.06 ±0.04 

0.65 ± 0.06 

0.13 ±0.02 


with an unknown exponent 7 which, as for the parameter /3 
in the previous section, presumably depends on the eccen¬ 
tricity of the cut-off ellipses. We therefore fit the extracted 
contours with the model 

r I J?, A, 7 ) = i? [1 -h A cos (uj cj))Y (24) 

and leave R free to vary while demanding that the param¬ 
eters A and 7 shall be equal for all curves. The results of 
fitting them simultaneously to the data are represented by 
the solid lines in Fig. 8 , with the corresponding best fit val¬ 
ues given in Tab. 2. 

We observe that the data is fitted well by the model, 
showing the expected distortion of radial scales, i.e. a 
stretching for those angles where the charactersitic scale is 
larger. Furthermore, the fitted values for R match well to 
the same isocontours of statistically isotropic fields, which 
were estimated from an ensemble of randomly orientated 
halos with fixed scales = 0.006 L and ay — a = 0.02 L. 
However, the parameter A is slightly off, which can be at¬ 
tributed to the scatter in the data and the fact that the fit 
is quite sensitive on both A and 7 . 


while keeping a^ = 0.0061/ fixed. As before, we intend to 
check how the ALCF is influenced by this anisotropy and 
especially how originally spherical lines of constant I are 
deformed. In order to suppress any additional variations of 
the amplitude of l{r±, ry), we need to hold the filling factor 
constant, so that 

fW = PH {(AS {(!>) = const. (21) 

=> PH{(t>) °^'^/S{(A ^ (22) 

where we resorted to the definition of / in Sec. 2.2 and used 
To = \/(^x <^y for D — 2. Consequently, when determining 
the scales from S, the respective orientation angles must be 
drawn from the distribution 1/S((f>)- This distribution is re¬ 
alised using the method of inverse transform sampling (De- 
vroye 1986). With that, we set up two ensembles of fields for 
the cases lo = 2 and w = 4, where we chose a = 0.02 L and 
Act = 0.01 L. After evaluating the line correlation function 
we extract two isocontours for each ui, convert the Carte¬ 
sian coordinates into 21 angular bins and average over 100 
realisations which gives the data points in Fig. 8 . 

Since the average characteristic scale ro = \Ja is 
modified by the angle-dependent factor -y/l -I- A cos {to cf) 
with A = It^afa, we may expect that any isocontour at 
scale R for a statistically isotropic field is distorted as 

i? [1-b A cos(tU((.)]'^ , (23) 


4 APPLICATION TO REDSHIFT-SPACE 
DISTORTIONS 

Anisotropies in the distribution of galaxies occur naturally 
whenever the redshifts and positions (angles) acquired in a 
survey are converted into physical distances. They may be 
of geometric or kinematic nature, denoted by the Alcock- 
Paczyhski effect and (dynamical) redshift-space distortions, 
respectively. In this section, we analyse whether these obser¬ 
vational effects are picked up by the ALCF and how these 
results compare to the ones from conventional two-point 
statistics. Since we primarily aim for a proof of concept, we 
will continue to use the two-dimensional halo and Zel’dovich 
mock fields introduced above. 


4.1 The Alcock-Paczyhski Effect 

The conversion of the redshift difference and angular sepa¬ 
ration of a pair of galaxies to a physical distance depends on 
two basic cosmological quantities: knowledge of the Hubble 
rate H{z) is required for the computation of the line-of-sight 
separation; and the angular diameter distance Da{z) is re¬ 
quired for the transverse separation. 



rj_ = Da{z) {1 + z) A9 . (26) 
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Consequently, creating a spatial map of the universe hinges 
on a cosmological model that has to be assumed a priori; 
if this model differs from the true cosmology, the deduced 
distances are wrong, causing distortions in the clustering 
signal. The first who proposed to use this effect as a probe 
for cosmological parameters were Alcock and Paczyriski (Al- 
cock & Paczyhski 1979) (AP). By considering an object 
whose intrinsic shape is known to be spherical, they deduced 
that measuring distortions from sphericity can be used to 
constrain the combination H{z) Da{z), and in turn, for in¬ 
stance, the equation of state parameter of dark energy. How¬ 
ever, in practice, there are no objects, which are sufficiently 
spherical or whose length scales are known to sufficient accu¬ 
racy, to perform the AP test on individually. Instead, one re¬ 
sorts to statistical standard rulers, that is, length scales that 
are statistically imprinted on large-scale structures, such as 
the one set by baryon accoustic oscillations (BAO); since 
BAO’s are primarily a linear phenomenon, their signal is 
not considerably affected by non-linear physics. Thus, by 
measuring the position of the BAO peak and comparing it 
to the theoretical prediction, it is possible to conduct an AP 
test (Blake & Glazebrook 2003). However, with the ever in¬ 
creasing volume coverage of galaxy surveys, an alternative 
has become feasible, in which a complete model of a given 
statistical measure, incorporating a possibly wrong choice of 
cosmology, is fitted to the data (see e.g. Ballinger, Peacock 
& Heavens 1996, Okumura et al. 2008 or Padmanabhan & 
White 2008). A similar method could be employed based on 
the ALCF. 

To quantify the impact of the AP effect on the line 
correlation function, we firstly parametrise the mismatch 
between true and assumed cosmology by two squashing fac¬ 
tors (following the notation of Ballinger, Peacock & Heavens 
1996): 


f Ha{z) 
^11 = ^ 


^ Da{z) 

DaAA ’ 


(27) 


where here and in the following the subscript ‘a’ indi¬ 
cates quantities derived from the assumed cosmology. From 
Eqs. (25) and (26) we see that assumed and true distance 
scales are related by the matrix 



Due to number conservation, the overdensity does not 
change in the assumed coordinate system, so that Sa{xa) = 
5{x), which implies for the two-point function 

dra) = CC = ra) ■ (29) 


Thus, any squashing of galaxy separations directly trans¬ 
lates to an equivalent squashing of the isocontours of To 
compute the analogue effect for the ALCF we need to be 
more careful because of the scale dependent cut-off and pref¬ 
actor. For that reason we transform la{r±,a,r\\A piecewise 
and start from the whitened and filtered density field in the 
assumed coordinate system, which we write as a convolution 


between e and the cut-off 0 in real space, 

ea * Oa,-r^{Xa) = Va j d^ fca 6*'““ 

^(fcaiT’a)^47r2 

= K|5r^ j Cfc. (30) 

6(S~^k,ra)^4:7r‘^ 

Here, we used that the Fourier modes scale inversely to spa¬ 
tial scales and employed the identity ta.fea = tfe, which is 
another direct consequence of number conservation. Shift¬ 
ing the 5-matrix from k to Va in the cut-off function causes 
an additional term that can be computed from the definition 
of 8 in Eq. (13), yielding 

e{S-^k, rC = e{k, 5-V,) + (/i - 4) (fci rl^) 

= eik, S-\a) + Spik, S-\a) , (31) 

SAk, r) = il- F^) {F-AI rl) , (32) 

where F is the ratio of /y and f± Plugging this result back 

into Eq. (30) gives 

Ca * Oa,r-a.iXa) = V J dke’^'^ ““ Cfc 

0(fc,<S~^T’a)+5j;’(fc,<S~^T’a)^47r2 

= e* 0F.s-i-r„(5~^a:a), (33) 

where we identified 14/151 with the true volume and let 0 f 
denote the altered cut-off that reverts to the original one 
for F = 1. Having thus determined the relation between the 
whitened and filtered density field in the assumed cosmology 
to the one in the true model, we still need to transform the 
prefactor. Using Eq. (28) again, we arrive at 


2 I 2 _ n —2 2 I n —2 2 _ / p — 2 ■ p — 2'\ 2 ■ p —2 2 

ra,\\ + ra,± = /|| r|| + r_L = (^/|| + f± ) ^|| + f± r 

= fz\^ [l+{F-^-l) A] , ( 34 ) 

where /r is the cosine of the angle to the line-of-sight, i.e. jj, — 
r\\/r. All in all, we get 

Lira) = Fi[l + {F-^ - 1) (5-V.)"] " If (5"V,) , 

( 35 ) 


where 


(5-Va)' 





2 

a,_L 


(36) 


Eq. (35), which is the main result of this section, relates 
the correlator for the assumed cosmology at scale Va to the 
modified (but true) one at the shifted position 5“^ r^. As 
in Eq. (29) for the two-point function, it encodes the impact 
of the AP effect on the ALCF and can either be used as a 
basis for fitting existing data to extract f± and /y or for 
simulating it. Since the effect is not immediately evident by 
inspecting Eq. (35), we compute the assumed line correlation 
function from a set of true fields with varying values for the 
squashing parameters. 


Note that 5p arises purely as a consequence of the first term 
in Eq. 13. Eq. 31 and all of the following expressions are hence 
equally valid for the original spherical cut-off. 
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Figure 9. Quadrupole-to-monopole ratio for the ALCF (upper 
panel) and two-point function (lower panel), derived from a set of 
50 Gaussian halo fields. The AP effect is simulated via Eq. (35) 
with /|| ranging between 1 and 1.5 while f± = l. Shaded areas 
mark l-cr uncertainties. 

To visualise the AP effect and enable comparison with 
the two-point function we decompose the data in the 
(fu, f|| )-plane into multipoles. As our density fields are two- 
dimensional, the n-th order multipole is calculated via 

/n(r) = — / f{r sin (j),r cos (j)) cos {n(l))d4>, (37) 

TT J-n 

where / is either ^ or I and the sign convention is cho¬ 
sen such that an enhancement of the transverse direction 
compared to the line-of-sight translates into a positive sig¬ 
nal for the quadrupole / 2 . Using Eq. (35), we determine 
the monopole and quadrupole for an ensemble of Gaussian 
halo and Zel’dovich density fields with varying /y, but fixed 
f±. The corresponding quadrupole-to-monopole ratios are 
plotted in the upper panels of Fig. 9 and 10. We notice 
that in both cases the AP effect introduces a clearly visi¬ 
ble anisotropy by means of an increasing quadrupole signal 
with rising /y. Interestingly, the quadrupole is positive for 
small scales, meaning that the ALCF is squashed. However, 
with increasing scale the quadrupole changes sign, turning 
the initial squashing into a stretching of the line-of-sight. Af¬ 
ter attaining a maximum, the amplitude of the quadrupole 
decreases, but remains negative throughout. In comparison, 
the quadrupole-to-monopole ratio of the two-point function, 
which is displayed in the lower panels of Fig. 9 and 10, 
is positive for all scales, reflecting the squashing expected 
from Eq. (29). Although being unintuitive, the behaviour 
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Figure 10. Same as Fig. 9, but based on 50 realisations of 
Zel’dovich density fields. 


of the ALCF can be understood in light of the discussion in 
Sec. 3.3. If the separations of galaxies are squashed in a given 
direction, the scale and orientation of the structures in the 
density field change. In the case considered here (squashing 
in the line-of-sight direction) a preference for structures to 
be aligned with the transverse direction arises. This results 
in the enhanced radial signal since the ALCF scales inversely 
with the number density of objects. 

While the ALCF is indeed capable of detecting the 
AP effect, the two-point function appears to be superior 
in telling apart different /y from each other for both sets 
of fields. Its quadrupole-to-monopole ratio increases with 
scale, which renders ^ sensitive to /y on all scales. On 
the other hand, the 1-a error regions begin to overlap at 
r ~ 0.02Gpc for the line correlation function, which 
implies that estimates based on ^ will be statistically more 
significant. This is a reasonable result owing to the num¬ 
ber of modifications in Eq. (35), which tend to influence 
the ALCF in converse directions. Assuming F > 1 in the 
following, we first note that the modified cut-off gives rise 
to an increase in the ALCF’s amplitude in the line-of-sight 
direction, as can be verified empirically. Conversely, the dis¬ 
tortion of the prefactor attains a minimum for ra,± = 0 and 
accordingly strengthens the transverse direction compared 
to the line-of-sight. The same is true for the modification 
of scales by the matrix 5“^, given that l{r) is a decreasing 
function of |r|. As a consequence, the net effect of the AP 
squashing on the ALCF is smaller than for the two-point 
function, where only one kind of distortion effect is present. 
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However, we would like to reiterate that the ALCF uses 
phase information, which is complementary to the ampli¬ 
tude information retained in the two-point correlator. Thus, 
both estimators can be used in conjunction in realistic ap¬ 
plications to tighten constraints. 

4.2 Kinematical Redshift-Space Distortions 

Another consequence of measuring redshifts instead of phys¬ 
ical distances are distortions due to the peculiar velocity of 
galaxies. Apart from the Hubble flow, they also contribute to 
the redshift and thus cause the galaxies to appear displaced 
along the line-of-sight. 

Since galaxies result from the structure formation pro¬ 
cess through gravity, their peculiar velocities are induced by 
gravity as well. On large scales, where structures have not 
yet fully collapsed, all matter tends to fall into the nearest 
overdensity. For that reason, a given galaxy that resides on 
the near-edge of a large cluster tends to move away from us, 
increasing its redshift. Consequently, its apparent position 
is closer to the center of the overdensity. A galaxy on the 
far-edge of the cluster would behave in the opposite way, so 
that it appears to be closer to us than it actually qis. We 
therefore observe that structures on those large scales look 
squashed in radial direction which is known as the Kaiser 
effect (Kaiser 1987). He showed that in the distant observer 
limit, the power of scales that are aligned with the line-of- 
sight is boosted by the factor I + ffT {fj, being the cosine 
to the line of sight), so that the redshift-space dark matter 
power spectrum is given by 

Ps{k) = {l + fti^fP{k). (38) 

The amplitude of the distortion depends on the logarithmic 
growth rate / = din D/d In a, which is related to the matter 
content in the universe and, according to Peebles (1980), well 
approximated by the power law / « Thus, measuring 
the Kaiser effect enables us to put constraints on the matter 
density, or, if combined with an independent measurement 
of Qm, on the theory of general relativity. 

At much smaller, non-linear scales, objects within al¬ 
ready collapsed structures acquire random virialized veloci¬ 
ties which smear the structure in the line-of-sight direction 
leading to a characteristic shape that seems to point at the 
observer. Somewhat misleadingly, these shapes were called 
Fingers-of-God (Jackson 1972). Lacking an analytical de¬ 
scription, several empirical models exist that try to mimic 
this small-scale effect by smearing the density field along 
the line-of-sight with the probability density function of ve¬ 
locities (Jackson 1972; Hawkins et al. 2003). The models 
mainly differ in the assumption of how these velocities are 
distributed; commonly, it is expected that they acquire a 
Maxwellian distribution, which gives rise to a multiplicative 
damping term for the power spectrum of the form 

DpoG(fc,/r) = ^1 H- ^■ (39) 

The pairwise velocity dispersion ap is approximately of order 
cTp ~ 400km/s (Hawkins et al. 2003). 

The goal for the remainder of this section is to es¬ 
timate how the ALCF is influenced by these kinematical 
redshift-space distortions. Hence, we take appropriate sets 
of mock fields that simulate the Kaiser or Fingers-of-God 
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Figure 11. Quadrupole-to-monopole ratio for the ALCF and 
two-point function calculated from a set of 50 Zel’dovich den¬ 
sity fields incorporating the Kaiser effect. The amplitude of the 
distortion is parametrised by /. 

effect and proceed as in the preceding section by comput¬ 
ing the quadrupole-to-monopole ratio. The Kaiser effect can 
be easily incorporated into the Zel’dovich density fields by 
boosting the Lagrangian displacement field in the line-of- 
sight direction by the factor 1 -I- /. This leads to distorted 
density fields that seem to be derived from redshift-space. In 
Fig. 11, we plot the quadrupole-to-monopole ratios for both 
the ALCF and the two-point function, estimated from these 
fields where we leave / as a free parameter which varies be¬ 
tween 0.25 and 0.75. While the two-point function displays 
the anticipated behaviour, i.e. a squashing in the line-of- 
sight direction that becomes more prominent for increasing 
/, the ALCF exhibits a more complex quadrupole signal. 
On small scales, the quadrupole is negative, indicating a 
stretching in the line-of-sight, before it changes sign at a 
crossover scale that seems to shift towards larger r for higher 
/. Thereafter, it remains positive, though the squashing is 
less pronounced than for the two-point function. Since the 
Kaiser effect shrinks the radial dimension of structures in 
the density field, the enhancement of the transverse direc¬ 
tion of the ALCF appears reasonable based on the discussion 
in Sec. 3.3.2. 

Since the Zel’dovich approximation cannot account for 
the Fingers-of-God effect, we create a simplified setup which 
mimics these small-scale distortions. Firstly, we draw ran¬ 
dom positions that are taken to be the centre points of ha¬ 
los. We assume that these halos have a spherical Gaussian 
density profile and accordingly draw particle positions in 
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Figure 12. ALCF quadrupole-to-monopole ratio for simple den¬ 
sity fields that mimic the Fingers-of-God effect. Each field con¬ 
sists of 100 spherical halos, each of which contains 1000 particles 
that are Gaussian distributed about the centre with standard de¬ 
viation 0.005 L. Their velocities are assigned via a Maxwellian 
distribution with velocity dispersion chosen such that the spread 
in position is cr = 0.015 L. Results are based on 50 realisations. 

each halo from a Gaussian distribution. Thereafter, we de¬ 
termine velocities from a Maxwellian distribution and dis¬ 
place all particles proportionally to their velocity along the 
line-of-sight. Computing the quadrupole-to-monop ole ratio 
gives the plot in Fig. 12, showing a distinctive negative sig¬ 
nal which confirms the stretching of structures in the line- 
of-sight. In this simplified case, the Fingers-of-God effect 
establishes a clear preference of the line-of-sight, with no 
structures being aligned with the transverse direction. How¬ 
ever, as we have seen in Sec. 3.3.1, the signal of the ALCF 
may change drastically if this is the case, so it would be 
interesting to see how the Fingers-of-God effect appears in 
a realistic N-body simulation that also contains filamentary 
structures in all directions. 


5 CONCLUSION 

In this paper, we presented an extension to the recently pro¬ 
posed line correlation measure (Obreschkow et al. 2013) that 
enables the identification of anisotropies in a given density 
field. The line correlation function is an estimator that re¬ 
lies purely on the phases of the density field and thus probes 
information not contained in ordinary two-point statistics. 

Starting from a generalisation of the originally defined 
line correlation function, we noticed that anisotropies tend 
to be smeared over the plane of radial and transverse sepa¬ 
rations. For that reason, we introduced a novel mode cut-off 
scheme, which plays an integral part in the definition of the 
anisotropic line correlation function. We showed that this 
new function produces better signal-to-noise ratios and, by 
using simple test fields as replacements for real observational 
data, we demonstrated that it extends the properties of the 
original line correlation function in a consistent manner. In 
particular, we observed that the altered cut-off method al¬ 
lows for a clear distinction between statistically isotropic 


and anisotropic fields and we showed that it is sensitive to 
angularly varying number densities and scales of structures. 

In the second part, but still aiming primarily for a 
proof of concept, we considered observational effects, like 
the Alock-Paczyhski, Kaiser and the Fingers-of-God effect. 
We derived the relation between the assumed and true line 
correlation functions when redshifts and angles of galaxies 
are transformed with a wrong prior cosmology and com¬ 
puted quadrupole-to-monopole ratios for all three effects, 
estimated from simple mock fields. We found that the line 
correlator is capable of detecting them, displaying different 
systematics than the two-point function. While its sensitiv¬ 
ity appears to be somewhat worse than that of the two-point 
function, this is still a promising result: since the information 
probed by both measures is independent from each other, 
it is conceivable that employing both statistics can lead to 
tighter constraints on cosmological parameters than, for in¬ 
stance, a combination of the two- and three-point function. 

To investigate these prospects, a more detailed study 
of redshift-space distortions is necessary. It is inevitable to 
employ N-body simulations to analyse both the Kaiser and 
Fingers-of-God effect in a realistic setup, and it would be 
interesting to see whether the crossover scale, which was 
observed in Fig. 11, is indeed a persistent feature. To un¬ 
derstand this behaviour, additional analytical studies are 
desirable as well. A possible ansatz for such an investiga¬ 
tion could be to reformulate the combination of bispectrum 
and power spectra that appears in the perturbative expres¬ 
sion of the line correlation function (see Eq. (3)) in terms 
of the perturbative kernels in redshift-space which were de¬ 
rived in Scoccimarro, Gouchman & Frieman (1999). Finally, 
it remains to be seen which parameters the line correlation 
function is able to constrain and to what accuracy. 

In addition to the independence of information probed 
by the two-point and line correlation functions, both mea¬ 
sures have different systematics, which allows the breaking 
of degeneracies in the parameter estimation. An interesting 
example is the logarithmic growth rate which describes the 
strenth of the Kaiser effect (see Sec. 4.2). Due to the fact 
that we are measuring galaxy positions, which are biased 
tracers of the underlying dark matter distribution, the two- 
point function is actually only sensitive to the combination 
/3 = f /b, where b is the linear bias parameter. On the other 
hand, the line correlation function is independent of linear 
bias and hence constrains / directly, meaning that a joint 
analysis might allow for the determination of / and b sep¬ 
arately. We leave these exciting new avenues for a future 
research project. 
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APPENDIX A: NOTES ON THE 
IMPLEMENTATION 

In this appendix, we briefly comment on our numerical im¬ 
plementation of the line correlation function and show con¬ 
vergence tests. 

We adopt a standard numerical discretisation scheme: 
the density field is contained in a finite cubic box of side- 
length L on which we impose periodic boundary conditions. 
The box is subdivided into an even number, N, of cells per 
side with the regular Cartesian grid spacing Aa; = L/N. 
Accordingly, in Fourier space the box has a sidelength 
2-17 N/L with spacing Ak = 2'r/L and Nyquist frequency 
fcNy = nN/L. Upon discretisation using this scheme, inte¬ 
grals in configuration and Fourier space become 

J d^xfix) ^/(a:), (Al) 

Jd^kU -r (A2) 

Applying these rules to the prescription of the line correla¬ 
tion function in Eq. (4) we get 

3D 

^(^■) = (^) " wodfei-fealr) 

|fel|.|fe2l, 

I+^2 I ^271 /r 

V ^hiSh2^—hi-k2 fA3l 

Implementing this expression is straightforward, but com¬ 
putationally expensive due to the nested summations over 
fei and fca. Hence, we use the equivalent of Eq. (4) in conhg- 
uration space, that is the spatially and rotationally averaged 
product of whitened density fields 



X eg{x + TZr) te{x — TZr))^ , (A4) 

where TZ denotes a rotation matrix and the mode cut-off is 
written as eg which stands for the convolution of e{x) with 
an appropriate hlter function Or- Discretising the equation 
above gives 






X ee{x + TZr) eg{x — 7?.r))^ 

(A5) 


3D 



y ^ FFT{(e,(a:) 



X eg{x + TZr) eg{x - TZr)).j ^}, 

(A6) 


where we replaced the summation over all x by the Fast- 
Fourier-Transform (FFT) algorithm. Since we can quickly 
transform back and forth between Fourier and conhguration 
space using the FFT algorithm, it is advantageous to com¬ 
pute the convolution as well as the shift in position in Fourier 
space. There, shifting the field by the vector r amounts to 
multiplication with the phase factor exp (ik ■ r), while min¬ 
imising rounding errors that might arise from the gridding 
when evaluating the held at a shifted position in conhgura¬ 
tion space. Hence, letting 9{k, r) either denote the original 
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or modified cut-off in Eq. (13), we have 

£ 0 ( 0 ;-I-r) = FFT“^ jcfc e*'” ’’6>(fe,r)| . (A7) 

Eqs. (A6) and (A7) form the basis of all numerical compu¬ 
tations in this article. 

Let us consider a test case that allows us to compare 
with analytic results and study the convergence of the im¬ 
plementation with increasing N. We create density fields 
consisting of Nh spherical Gaussians with equal standard 
deviation a. In this case (see Eqs. (9) and (10)), 

5(fcl,fc2,fc3) ^ / V 

VP(|fel|)P(|fc2|)P(|fc3|) V ’ 

such that for D = 2 and, after integrating out one azimuthal 
angle, Eq. (3) gives 

0 0 

Mcut 

/ - 7 ^=M\ki-k 2 \r), (A9) 

4 V 1 - 

where fj, = cos<(fci,fc 2 ) and its upper limit is given 
by Mcut = min{l,max{-l, [(27r/r)^ - kf - A:i]/[ 2 A:ifc 2 ]}}. 
Eq. (A9) can be integrated easily, for instance using a 
Monte-Carlo approach, to yield the analytic answer. How¬ 
ever, to compare with the numerical result we need to take 
care of a subtlety which arises in the whitening process: 
whenever the amplitude of Fourier coefficients drops below 
a certain threshold value, their phase factors are not re¬ 
solved properly any more; in order to avoid the introduction 
of artifacts, we need to cut out these coefficients®. We can 
apply the same rule as to the zero mode and set all modes 
with |(5fc| < 10“^ to zero. For our test fields there exists an 
isotropic threshold scale fcth so that all modes k > fcth have 
unresolved phase factors. Hence, this effect can be dealt with 
in Eq. (A9) by replacing the integration limits for k\ and fc 2 
by min{fcth, 27r/r} and similarly in /Tcut. 

Another way to deal with unresolved phase factors is the 
superposition of the density field with a Gaussian random 
field (GRF) of small amplitude to ensures that Fourier co¬ 
efficients with nearly vanishing amplitudes acquire random 
phases. This procedure guarantees that no artificial infor¬ 
mation compromises the signal: since the Gaussian random 
field is uncorrelated to the Gaussian halos and assuming 
that it has zero mean, we only obtain a contribution to the 
overall power spectrum, 

P{k) = Paik) + PGRF{k). (AlO) 

Doing so complicates Eqs. (A8) and (A9), but for known 
PcRpik), it is still possible to find the analytic answer via 
Monte-Carlo integration. 

The upper panel of Fig. Al shows a comparison be¬ 
tween numerical and analytic results. We see that to the 
right of the vertical line, which indicates the threshold scale 
Tth = 27r/fcth, the data points approach the dashed line, 



r[L] 


Figure Al. Top panel: Comparison between analytics and nu¬ 
merics for a density field comprising Nr = 100 spherical Gaussian 
halos (cr = 0.024L). The dashed line is the result from Eq. (A9) 
while the blue and green solid lines either account for a fixed 
threshold scale or the superposition with a Gaussian random field 
with Pgrf = const.; the corresponding data points stem from the 
purely numerical implementation averaged over 50 realisations. 
The vertical line indicates the threshold scale rtj, = 2f Bot¬ 
tom panel: Dependency on the number of grid points for a single 
realisation of the same density fields. 

which is the analytic result based on Eq. (A9). Below rth, 
the effect of cutting out unresolved phase factors becomes 
important and produces the oscillatory behaviour, which can 
be exactly reproduced for both methods if we account for 
the corresponding modifications in Eq. (A9). Thus, this test 
case serves as a convincing consistency check between our 
numerical scheme and the perturbative expansion of the line 
correlation function. 

Finally, the lower panel of Fig. Al displays a conver¬ 
gence test for a varying number of grid points using a single 
realisation of the test fields above. The plot clearly shows 
that our method converges quickly, indicating only small 
discrepencies for N as low as N = 128. All of our computa¬ 
tions were carried out with N = 512 grid points. 


® This feature is due to the peculiar nature of the test fields 
considered here. Realistic density fields are unlikely to exhibit 
unresolved phase factors. 
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